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Topological defects play a pivotal role in the physics of liquid crystals and represent one of the most 
prominent and well studied aspects of mesophases. While in two-dimensional nematics, disclinations 
are traditionally treated as point-like objects, recent experimental studies on active nematics have 
suggested that half-strength disclinations might in fact possess a polar structure. In this article, 
we provide a precise definition of polarity for half-strength nematic disclinations, we introduce a 
simple and robust method to calculate this quantity from experimental and numerical data and we 
investigate how the orientational properties of half-strength disclinations affect their relaxational 
dynamics. 


I. INTRODUCTION 

Disclinations are singularities in nematic liquid crystals 
where the average molecular orientation is undefined [1- 
3]. These defects are ubiquitous in nematic samples and 
represent their most prominent and visible feature [4-10]. 
In fact, the existence itself of a nematic phase, where the 
molecules are oriented along a common direction, but 
have zero macroscopic polarization, was established only 
once disclinations were correctly identified [11]. 

In two-dimensional nematic liquid crystals, disclina¬ 
tions are conventionally treated as point-like objects in¬ 
teracting with each other through elastic forces formally 
analogous to the Coulomb force in electrostatics [1-3, 12- 
lb]. Recently, however, combined experimental and theo¬ 
retical efforts toward understanding the mechanics of ac¬ 
tive nematics (i.e. nematic liquid crystals obtained from 
self- or mutually-propelled rod-like macromolecules, typi¬ 
cally of biological origin), have suggested that some type 
of nematic disclinations have in fact a polar structure 
and can reorient each other via elastic and hydrodynamic 
torques [17-20]. Two spectacular examples of these phe¬ 
nomena have been recently reported in two-dimensional 
active nematic suspensions of microtubule bundles and 
kinesin. Keber et al. [17] constructed an active ne¬ 
matic vesicle by encapsulating microtubules, kinesin mo¬ 
tors clusters and polyethylene glycol (PEG) within a 
lipid vesicle. A depletion mechanism, due to the PEG, 
drives the microtubules to the inner leaflet of the vesicle, 
giving rise to a dense two-dimensional nematic cortex. 
As a consequence of the spherical confinement, such a 
two-dimensional nematic is forced to contain four +1/2 
disclinations [21-23], which, due to the local hydrody¬ 
namic flow fueled by the active stresses, travel at con¬ 
stant speed toward the “head” of their comet-like struc¬ 
ture [17]. More recently, DeCamp et al. [18] demon¬ 
strated that defects in two-dimensional active nemat¬ 
ics can themselves form a nematic phase in which the 
head-to-tail directionality of +1/2 disclinations propa¬ 
gates over distances several order of magnitude larger 
than the length of a single microtubules [18]. 

In spite of this convincing experimental evidence, a 
precise definition of defect orientation is, however, still 
lacking. Intuitively, the existence of an orientation per 


se does not depend on whether the system is passive or 
active, but it is only determined by the defects local ge¬ 
ometry. Here we consider the simplest possible setting, 
consisting of a two-dimensional passive nematic liquid 
crystal subject to a purely relational dynamics, and we 
address the following fundamental questions: what is be¬ 
hind the comet-like appearance of +1/2 disclinations? 
Can an orientation be defined for —1/2 disclinations as 
well? How can the orientation of a defect be determined 
from numerical and experimental data and, perhaps more 
importantly, does such an orientation affect the mechan¬ 
ics of half-strength disclinations in any way? 

We start by providing a rigorous definition of defect 
orientation and introduce a simple method to calcu¬ 
late the orientation of half-strength disclinations from 
pixelated data. Next, using numerical simulations, we 
demonstrate that pairs of like-sign disclinations exert 
elastic torques that tends to anti-align them as they re¬ 
pel, while oppositely charged disclinations have negligible 
orientational interaction. For both cases, however, the 
coupling between translational and rotational dynamics 
is very strong and leads to a variety of novel annihila¬ 
tion/repulsion trajectories. Finally, we analytically ad¬ 
dress the ideal case of two like-sign disclinations on an 
infinite plane and show that these interact through an 
elastic torque independent on their distance and remi¬ 
niscent of that between dislocations in two-dimensional 
solids [24, 25]. As a consequence, pairs of like-sign discli¬ 
nations on an infinite plane anti-align exponentially with 
time. 


II. RESULTS 

A. Polarity of half-strength nematic disclinations 

Let us consider a nematic liquid crystal in two dimen¬ 
sions and let n = (cos 0, sin 0) be the nematic direc¬ 
tor representing the average molecular orientation. In 
the presence of disclinations, the director rotates by a 
multiple of i r in one loop around the defect core, thus: 

d6 = 2irk , where fc = ±1/2, ±1 ... is the turning num¬ 
ber of “strength” of the disclination and the contour in¬ 
tegral is calculated on any closed path enclosing the core. 
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In one elastic constant approximation, a well-known de¬ 
fective solution that minimizes the Frank energy: 

e f = l -K jdA\ve\ 2 , (1) 

with K an elastic constant, is given by [1]: 

0 = k(t> + 0 o , ( 2 ) 

where (p = arctan (y/x) is the usual polar angle. The con¬ 
stant $o, which is often omitted, describes a global rota¬ 
tion of the director about the z—axis passing through 
the defect core. Now, as it is evident from Fig. 1, 
these energy minimizing defective configurations have 
2 1 1 — k\— fold rotational symmetry. This can be made ex¬ 
plicit by introducing a new constant angle ip, such that 
#o = (1 — k)ip. Eq. (2) becomes then: 

0 = k(<p - ip) ± 'ip • (3) 

For k = 1/2, 0 = (</> + ?/)/2 and it takes a full 2tt rotation 
in order to transform n into itself: i.e. n{^ip) = n(ip + 
2m7r), with m an integer. While for k = —1/2, 6 = 
— ((/)— 3 , 0)/2 and n is manifestly invariant under rotations 
by multiples of 2tt/3: i.e. n(ip) = n(/0 + 2ra7r/3). The 
vector: 


p = (cost/, sin'll;) , (4) 

with | ip | < 7r/[2(l — k)] defines then the polarity of discli- 
nations of strength k = ±1/2 (see Fig. 1). By construc¬ 
tion, p is defined up to rotations by 7r/(l — k). 

We now turn to the question of calculating the defect 
orientation from an arbitrary defective configuration of 
the nematic director. Both in experiments and numerical 
simulations, local orientations are generally available in 
the from of a discrete set of data points on a grid, say 
{ 0i }, with i the label of a generic grid point. Disclinations 
can then be detected by measuring the turning number 
of each plaquette V in the grid [26]. Then, the simplest 
method to track the orientation consists of constructing 
a sequence of data pairs (0^, (pi), with i G P and (pi mea¬ 
sured with respect to the center of the plaquette, and us¬ 
ing directly Eq. (3) to find ip through a linear regression. 
This method, which was adopted in Ref. [18], is however 
often unsatisfactory. Eq. (3) corresponds indeed to the 
ideal, energy minimizing, configuration of a disclination, 
while in practice the nematic director will always appear 
distorted. As a consequence, the data points (0$, (pi) will 
never align along a straight line, but rather form an un¬ 
dulated curve (see the supplementary information of Ref. 
[18]), thus introducing the risk of uncontrolled systematic 
errors. 

In order to overcome this difficulty, we introduce an 
alternative method directly inspired to the mechanics of 
disclinations in active nematics [17, 27]. In Ref. [27] 
it was demonstrated that active ±1/2 disclinations self- 
propel in the direction of their symmetry axis by virtue of 
the spontaneous flow powered by the active stresses. The 



FIG. 1. Examples of +1/2 (a) and —1/2 (b) disclinations with 
generic orientation p = (cos 'ip, sin?/;), with vp calculated from 
Eq. (5). For —1/2 disclinations, p is related to the polarity 
p' — (cos?//, sin ip') of the dual +1/2 disclination obtained 
from the mirror-reflection 0 — > —0, as ^p = —ip'/3. 


body force driving the active flow is: / a oc V Q = p/ (2 r) 
with Qij = S(riirij — Sij/2) the nematic tensor [1], S 
the order parameter and r the distance from the de¬ 
fect core. The polarity of +1/2 disclinations can then 
be simply calculated from the divergence of the nematic 
tensor: p = V- Q/|V-Q|. In order to extend this 
method to —1/2 disclinations, we notice that +1/2 and 
— 1/2 disclinations can be transformed into one another 
by a mirror reflection of the director: 0 -+ —0. The 
polarity p' = (cos ip', sin?/') of the dual +1/2 disclina¬ 
tion, obtained by mirror-reflecting a —1/2 disclination, 
is straightforwardly related to that of the original —1/2 
disclination as ip m —'ip'/ 3 (Fig. 1). This allows us to cal¬ 
culate the angle pj for both positive and negative discli¬ 
nations: 


k (sgn(fc) d x Q X y dyQ xx ) 

- -- arctan ——- y - y —- 

1 -k L {OxQxx ± Sgn(fc) OyQxy) 


( 5 ) 


where (• • •) denotes an average along the shortest avail¬ 
able loop enclosing the core and we used the fact that 
Qxy -+ —Qxy under mirror reflections, hence the sign 
function sgn(fc). Eq. (5) provides a simple and robust 
way to calculate the orientation of a defect that is suit¬ 
able to be implemented on structured and unstructured 
grids upon approximating the derivatives by finite dif¬ 
ferences and using a plaquette to calculate the average. 
Supplementary Movie SI shows an example of this algo¬ 
rithm in a system undergoing a fast coarsening from a 
random configuration. 

Evidently, the auxiliary vector p', used in the deriva¬ 
tion of Eq. (5), can be also employed to calculate the 
relative orientation of two —1/2 disclinations. It is worth 
to stress, however, that p' is not uniquely defined as its 
direction depends on the choice of the axis used for the 
mirror transformation (the x —axis in this case). Choos¬ 
ing a different axis, say x r = (coscqsum), yields a dif¬ 
ferent vector p' with ?// = —3 ('ip — a). This ambiguity, 
however, has no effect on the relative orientation of an 
arbitrary number of —1/2 disclinations. Furthermore, as 
the Frank free energy is invariant with respect to mir- 
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FIG. 2. (a) The energy difference A E = FTdG(A'0) — 

F'LdG (7r), with F/LdG given in Eq. (6), of a pair of +1/2 
disclinations as a function of their relative tilt At/ for var¬ 
ious distances d of the cores. Inset: the same data rescaled 
by L/d and plotted versus (At/ — 7r) 2 , showing a dependence 
of the form A E ~ KL/d(A^p — 7r) 2 . (b) Two +1/2 disclina¬ 
tions with At/ = 7t/2. The rotation of the director within the 
transition region (shaded) between the cores, determines the 
angular dependence of the elastic energy. 



ror reflection, the energy of an arbitrary distribution of 
— 1/2 disclinations is equal to that of a distribution of 
+1/2 disclinations having the same positions and rela¬ 
tive orientations. The consequence of this invariance will 
be analyzed in the following sections. 


FIG. 3. The product p\ • p 2 as a function of time for a pair 
of +1/2 (a) disclinations and the corresponding trajectories 

(b) in the xy— plane. The data are obtained from a numerical 
integration of Eq. (7) with e/L — 5 x 10 -3 , initial distance 
d/L — 0.05 and various initial tilt values. The product p' x p 2 

(c) and the corresponding trajectory (d), for a ±1/2 pair. 


B. Orientational energy 


Now that a notion of orientation has been precisely 
defined and a method to calculate it has been introduced, 
we can examine how this affects the mechanics of half¬ 
strength disclinations. With this goal, we introduce the 
Landau-De Gennes energy [1, 28]: 


^LdG — 


\ K ! 


dA 


|VQ| 2 


■ tiQ 2 (trQ 2 - 1) 


(6) 


with e a constant with dimensions of length and propor¬ 
tional to the defect core radius. Eq. (6) is better suited 
for numerical applications than the Frank free energy as, 
unlike the nematic director n, the Q tensor is defined 
also at the defects (where the order parameter S drops 
to zero). In order to quantify the energetic cost of de¬ 
fect orientation, we consider two +1/2 disclinations in a 
square L x L domain, positioned at r% = (—d/ 2,0) and 
7*2 = (d/ 2,0) and having t/i = tt and varying 7 / 2 - Intu¬ 
itively, we might expect the antiparallel configuration for 
which At/ = t/i — 7/2 = tt to be energetically favorable. 
Fig. 2a shows the difference A E = E^q(A'iP) — E^q(tt) 
as a function of the angular displacement At/ and for var¬ 
ious distances d. The data show a clear dependence of 
the form A E ~ KL/d (At/ — 7 r) 2 (inset). 

The origin of this behavior is not difficult to under¬ 
stand and provides important insights in the orienta¬ 


tional mechanics of nematic disclinations. A pair of +1/2 
disclinations consists of three regions: the regions sur¬ 
rounding the cores where 0 = arctan[(T/ — yi)/(x — Xi)\/2 
with i = 1, 2, and an intermediate transition region 
(Fig. 2b). Within the transition region, the director 
is forced to rotate by an amount proportional to tt — At/ 
in order for 6 to match the orientation of the core re¬ 
gions. If the defects are separated by a distance d, the 
area of the transition region is roughly A b ^ Ld while 
the energy density scales like e^/K ~ (At/ — it) 2 /d 2 . 
Thus the total energy introduced by the transition re¬ 
gion is e\yA\y ~ (KL/d)( At/ — 7 r) 2 . In the general case, 
the energy of a pair of like-sign defects will depend on 
the detailed geometry of the transition region (i.e. how 
smoothly or abruptly the director rotates while interpo¬ 
lating between the cores), but the energy scaling will re¬ 
main unchanged. 


C. Relaxational dynamics of disclination pairs 

If the disclinations are now left free to relax, we ex¬ 
pect them to rotate and progressively reach the energy¬ 
minimizing antiparallel configuration, thus effectively ex¬ 
erting on each other an elastic torque T = —dE/dip. In 
order to test this scenario, we have numerically integrated 
the equation governing the relaxational dynamics of the 
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FIG. 4. The temporal evolution of a pair of +1/2 (a-c) and 
±1/2 (d-f) disclinations obtained from a numerical integra¬ 
tion of Eq. (7) with e/L = 5x 10 —3 , initial distance d/L = 1/3 
and initial tilt A if; — 0. 


nematic phase: 


dQij _ 1 SE Ld Q _ AT T 1 -S 2 

dt ~ 7 dQij ~ 7 r ViJ * 2 


( 7 ) 


with 7 the rotational viscosity [1]. Eq. (7) has been 
integrated using finite differences on a 256 x 256 grid 
and Neumann boundary conditions: i.e. d±Qij = 0 
at (x,y) = (±L/ 2 t ±L/ 2 ), with d± the normal deriva¬ 
tive (see Supplementary Movies S2-4). Choosing a fixed 
orientation at the boundary or periodic boundary con¬ 
ditions would restrict the type of defect pair we could 
consider, as a consequence of the Poincare-Hopf theo¬ 
rem [29]. Neumann boundary conditions reflects the ex¬ 
perimental scenario of weak or no anchoring. Within 
Landau-de Gennes theory, surface anchoring is often de¬ 
scribed via the Nobili-Durand energy density [30] given 
by: eND = W tr(Q — Q o) 2 , where W is a constant and 
Qo embodies the preferred orientation at the boundary 
(typically tangential or orthogonal). Neumann boundary 
condition corresponds to the case in which W ~ 0 and 
the director at the boundary is free to rotate. 

Fig. 3a shows the dot product p\ • P 2 as a function of 
time for two + 1/2 disclinations with various initial tilt. 
Regardless of the initial orientation, the defects relax to¬ 
ward the antiparallel configuration pi = —p 2 . The re¬ 
laxation rate is not constant and depends on the distance 
from the boundary, which becomes progressively shorter 
as the defects repel. For short times, while the defects 
are still far from the boundary, the angular displacement 
decays exponentially as one would expect from energetic 
considerations. 

The angular dynamics is strongly coupled with the 
translational dynamics of the core. The latter has been 
described by Kawasaki [31] and Denniston [32] and re¬ 
lies on the decomposition 6 = 6 1 + 0 ext , where 6\ = 
ki arctan [{y — yi)/(x — x-\)] describes the behavior near 
the core at rq = (xi,yi) and # ex t describes the depar¬ 


ture from the optimal defective configuration and plays 
the role of an external field. If the core radius is suffi¬ 
ciently small, the dynamics of the core is governed by the 
equation: 


dri 

c-^ = -27 ThKV ± 9 ext ( ri ) , ( 8 ) 

with Vi = (— d x ) [32]. The quantity £, is an effective 
drag coefficient proportional to the rotational viscosity 
7 in Eq. (7) and generally dependent on the core ve¬ 
locity: £ « 7ry/c 2 log(3.6/Er), with Er = ^a\dr\/dt\/K 
is the core Eriksen number. In first approximation 
log(3.6/Er) « 1 as a defect typically moves by a few 
core radii within the nematic relaxational time scale 
r = 7 a 2 /iF, hence \dr\/dt\ ~ a/r. 

In the presence of a second defect located at r 2 = 
(^ 2 , 2 / 2 ), approximating 6> ext = k 2 aictaii[(y-y 2 )/(x-x 2 )] 
and using (8), yields V_L0 ext (ri) = -k 2 r 12 /\ri 2 1 2 , with 
7*12 = rq — r 2 . Thus pairs of disclinations attract or re¬ 
pel each other as over-damped charged particles subject 
to the two-dimensional Coulomb force F\ 2 = —F 2 \ = 
2i:kik 2 Kri 2 /\ri 2 \ 2 . Fig. 3b shows the trajectories of the 
two defects in the central region of the square domain. 
The additional distortion due to the departure from the 
anti-parallel configuration introduces a transverse com¬ 
ponent in Vj_ $ext •> causing the defect to follow a curved 
trajectory (Fig. 4a-c). Everything we discussed for pairs 
of +1/2 disclinations holds for pairs of —1/2 disclinations 
as well since the energy (6) is invariant under the mirror- 
reflection (i.e. Q xy -+ — Qxy) that transforms +1/2 into 
— 1/2 defects. 

The case of a ±1/2 pair appears, on the other hand, 
very different as demonstrated by the dynamics of the 
product p' x • p 2 shown Fig. 3c. In this case, the defects 
maintain the same orientation during the entire annihi¬ 
lation process, while the trajectories of the cores reveal 
the same prominent angular dependence observed for the 
+1/2 pair (Fig. 3d and Fig. 4d-f). 


D. Orientational dynamics of like-sign disclination 
on an infinite plane 


From the previous section it should be clear that a 
pure particle description of the rotational dynamics of 
± 1/2 disclinations is not possible because the energy of 
a given configuration, hence the torque, depends on the 
structure of the transition region between the defects. A 
special situation, where this limitation does not occur, 
is represented by the case of two like-sign disclinations, 
located on an infinite plane at 77 and r 2 and oriented 
along the directions pi and p 2 . This problem, which was 
preliminarily discussed in Ref. [17], can be addressed 
by introducing an oppositely charged “image defect” at 
r* sa Vi + Rpi, with i = 1,2 and R an arbitrary distance 
(Fig. 5). The energy of this auxiliary configuration is 
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FIG. 5. Illustration of the method used for the derivation of 
Eq. (10). In order to introduce the energetic contribution 
due to the relative orientation of two like-sign defects, a pair 
of +1/2 disclinations (a) is replaced by a pair of ±1/2 dipoles 
(b), such that r* = ri + Rpi (i — 1 , 2), with ri (r*) the 
position of the positive (negative) defect and R an arbitrary 
distance. The energy of the original configuration can be 
recovered by taking R —> oo. 


given by: 


1,4 

E = -2t tK E kikj log 

i<j 



( 9 ) 


with a the defect core radius. The energy of the original 
configuration can be then recovered by taking R oo 
while preserving the orientations. The yields, after sim¬ 
ple algebraic manipulations: 


E = E/seif — 27rk 2 K 


log 


n - r 2 


2 lo g(! “ Pi • P2) 


( 10 ) 

where £7 se i f rsj log (R/a) is a position and orientation 
independent self-energy. The orientational contribution 
is then minimal when pi = —P2, while it diverges for 
pi = P 2 - More remarkably, the second term in Eq. (10) 
does not depend on the distance between the defects. 
This latter property, which is reminiscent of the behavior 
of dislocations in two-dimensional solids [24, 25] is due to 
the fact that, in this construction, the transition region 
between the two disclinations occupies the entire plane. 
Thus, while the positional interaction still depends on 
the defect separation, the orientational interaction is de¬ 
localized over the entire plane. The torque associated 
with Eq. (10) can be readily calculated: 


T12 = -T 21 = -7 Tk 2 K (11) 

1 - Pi • P2 

where pf- = (—p y: p x )- Consistently with Eqs. (7) and 
(8), we postulate the orientational dynamics resulting 
from the torque (11) to be purely overdamped, so that: 


t (hjH 
^ dt 


= 7 xk 


■K cot 

i<j 




( 12 ) 


with ( r an effective rotational drag coefficient. Using Eq. 
(12) one straightforwardly finds: 


cos A 2 p(t) + 1 = [cos A^O) + 1] e * , (13) 

with r = £ r /(27r Kk 2 ). Thus the defects exponentially 
relax toward the anti-parallel configuration, consistently 
with the early times dynamics observed in our simula¬ 
tions. The construction out-lined above cannot be ex¬ 
tended to the case of a pair of ±1/2 disclinations as a 
consequence of the non-uniqueness of the p' vector. 


III. DISCUSSION AND CONCLUSION 


Although discovered in the context of active nematics 
[17, 18], the polar structure of half-strength disclinations 
is a general property of nematic defects in both passive 
and active systems. In this article we demonstrated that 
a notion of polarity can be introduced for both +1/2 and 
—1/2 disclinations via the the vector p defined in Eq. (4) . 
Due to the discrete rotational symmetry of half-strength 
defects, p is defined up to rotations by 7r/(l — fc), with k 
the defect turning number. Thus, for +1/2 disclinations, 
p spans the entire unit circle, while in the case of —1/2 
disclinations it is defined up to rotations by 27 t/ 3 (Fig. 
1). As in the case of dislocations in solids, the elastic en¬ 
ergy of a pair of like-sign disclinations, depends on their 
relative orientation, with the antiparallel configuration 
representing the lowest energy alignment (Fig. 2). As a 
consequence, half-strength disclinations effectively exert 
torques on each other and, if left free to move, exponen¬ 
tially relax toward the antiparallel configuration (Figs. 3 
and 4). 

While in this article we have laid down the first fun¬ 
damental concepts of the orientational mechanics of ne¬ 
matic disclinations, much is still to be understood. Both 
in passive and active nematics, the strong distortion as¬ 
sociated with defects fuels hydrodynamic flow, which 
in turn causes a reorientation of the nematic director 
[14, 27, 33, 34]. The structure of such a backflow , is deter¬ 
mined by the geometry of the defect, thus by its polarity. 
This suggests that, in addition to the elastic torques de¬ 
scribed here, half-strength disclinations can additionally 
exert hydrodynamic torques. These are expected to lead 
to a richer and more complex type of orientational inter¬ 
actions, especially in the case of active nematics where 
backflow effects are more pronounced. 
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